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We study stationary electron flow through a three-terminal quantum ring and describe eflfects due 
to deflection of electron trajectories by classical magnetic forces. We demonstrate that generally at 
high magnetic fleld (B) the current is guided by magnetic forces to follow a classical path which 
for B > Q leads via the left arm of the ring to the left output terminal. The transport to the left 
output terminal is blocked for narrow windows of magnetic field for which the interference within the 
ring leads to formation of wave functions that are only weakly coupled to the output channel wave 
functions. These interference conditions are accompanied by injection of the current to the right 
arm of the ring and by appearance of sharp peaks of the transfer probability to the right output 
terminal. We find that these peaks at high magnetic fleld are attenuated by thermal widening 
of the transport window. We also demonstrate that the interference conditions that lead to their 
appearance vanish when an elastic scattering within the ring is present. The clear effect of magnetic 
forces on the transfer probabilities disappears along with Aharonov-Bohm oscillations in a chaotic 
transport regime that is found for rings whose width is larger than the width of the channels. 

PACS numbers: 73.63.-b, 73.63.Nm, 73.63.Kv 



I. INTRODUCTION 

Phase-coherent electron transport in mesoscopioi"— 
and nanoscale^^— rings results in appearance of 
Aharonov-Bohni^ conductance oscillations in external 
magnetic field. These conductance oscillations arc 
extensively studied in the context of scanning gate 
spectroscopy)^ spin-orbit coupling for both electronsiiii^ 
and holeSfi^ Aharonov-Bohm interferometrj*i^ includ- 
ing electron self-interferenceyi^ violation of Onsager 
symmetryfi^ and magnetic forces J^"— 

The deflection of electron trajectories by magnetic 
forces in two-terminal quantum rings was previously 
studied by time-dependent wave packet simulationSiii 
which indicated that in presence of external perpendicu- 
lar magnetic field the electron packet is preferentially in- 
jected into one of the arms of the ring which reduces the 
Aharonov-Bohm interference of electron waves meeting 
near the exit to the output lead. A time-dependent sim- 
ulation was also used to describe the transport through 
a three-terminal quantum ring^i^ which demonstrated 
that the Lorentz force - besides the reduction of the 
Aharonov-Bohm oscillations at high field - results in a 
distinct imbalance of the wave packet transfer probabili- 
ties to the two output leads. Such an imbalance of con- 
ductance of two output leads was indeed found in a recent 
experiment 1^ 

A three-terminal quantum ring is a basic elemental for 
construction of ring arrays that are proposed for imple- 
mentation of quantum logic operations22i2^ using spin- 
orbit interactions. In these structures^ir— the direction 
of the charge current is determined by the electron spin 
orientation. Magnetic forces^ can provide a mean of ex- 
ternal control of the current flow. 

The time-dependent simulations as previously per- 



formed for three-terminal ring a^^i^^ are based on a rela- 
tively straightforward procedure that indicates in a clear 
way the electron trajectories across the nanostructure. 
The charge transfer is a time-dependent process only 
in selected experiments, cf. the single-electron injec- 
tion into the quantum ringi^ realized according to the 
single-electron pump technique based on the Coulomb 
blockadci^ The standard experiments measure the cur- 
rent due to the stationary electron flow at the Fermi level 
which is therefore of a more basic interest than the wave 
packet dynamics. With the time-dependent approach 
one can in principle approach the monoenergetic time- 
independent limit increasing the spatial spread of the 
wave packet in the initial condition but the latter is lim- 
ited by necessarily flnite size of the computational box. 

The purpose of the present paper is to describe the ef- 
fect of magnetic forces on electron transport through a 
three-terminal ring in Hamiltonian eigcnstates. We find 
that at high magnetic field the electron flow follows the 
path determined by the Lorentz force - one of the arms of 
the ring is selected by the current which leaves the ring 
to the nearest output channel. However, exceptions to 
this rule are found for resonant interference conditions 
that block the transport to the output channel that is 
preferred by magnetic forces. This blockade is accompa- 
nied by anomalous (nonclassical) injection of the current 
to the ring and by appearance of peaks of the transfer 
probability to the other output channel. We study the 
thermal stability of this anomalous current injection, the 
influence of the elastic scatterers for the resonant inter- 
ference and effects of magnetic forces in chaotic trans- 
port regime. We also study oscillations of the current 
circulation which turn out to be more thermally stable 
than the oscillations of the transfer probabilities. Ori- 
entation of the currents circulating inside the ring de- 
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termines the sign of the magnetic dipolc moment that 
they generate. The magnetization oscihations due to the 
Aharonov-Bohm effect were so far measured for mcso- 
scopic open quantum rings2^ and for large ensembles of 
closed nanoringsi^ 



II. THEORY 

A. Model system 

The geometry of the studied system is depicted in Fig. 
[TJ The inner and outer radii of the ring are 88 nm and 
154 nm, respectively. The channels are assumed 68 nm 
wide. We treat the straight channel connected to the ring 
from below as the input terminal. The contacts to the 
input and the output channels are spaced by 120° angles. 
The output channels are bent twice under the angle of 
30° to acquire vertical orientation at the end of the com- 
putational box which allows for a uniform treatment of 
incoming and outgoing wave functions and currents (see 
below). 

We adopt a two-dimensional model and assume that 
the magnetic field is oriented perpendicular to the plane 
of confinement. We consider the electron Hamiltonian in 
form 



H 



{-ifiV + eA{r)Y /2m* + V{x, y) 



(1) 



where V{x, y) is the confinement potential - assumed 
zero within the channels (white area in Fig. [T|) and 
Vq = 200 mcV in the outside (the grey area in Fig. [1]). 
The potential offset Vq corresponds to channels made of 
GaAs embedded in an Alo.45Gao.55As matrix. In Ec^. ([T]) 
— e is the electron charge (e > 0) and m* = 0.067too is 
the GaAs electron band effective mass. 
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FIG. 1: Schematic drawing of the three-terminal ring. The 
confinement potential is zero inside the channels and 200 meV 
in the outside (grey area). The channel width is 68 nm, the 
inner and outer radii of the ring are 88 nm and 154 nm, re- 
spectively. The channel connected to the ring from below is 
the input lead. 



yields 



2m 



^ (p + eA)'^'^,, 



4*„ 



- . (2) 



- "^(Xfj^^y^), Cy = exp [-i|Aa;Ay], and 



where ^ 

Cx — exp [— i^AxAa;] . We apply the Lorentz gauge 
A — [Ax, Ay, O) = (O, Bx, O), for which the mesh Hamil- 
tonian eigenequation reads 



B. Hamiltonian discretization 

For the description of the stationary charge transport 
through the system we need to determine the Hamil- 
tonian ([T]) eigenfunctions for the electron coming from 
the input channel. We employ the finite difference ap- 
proach with a square computational box of side length 
482 nm (see Fig. [T|) on a grid of 241 x 241 points with 
mesh spacings Ax = Ay = 2 nm. The results pre- 
sented below are unaffected when one enlarges the com- 
putational box to cover a larger part of the input and 
the output channels. We use the Wilson^i type of dis- 
cretization of the kinetic energy operator in a version 
adapted by Governale and Ungarelli^^ for semiconduc- 
tor nanostructures. The discretization is consistent with 
the original Hamiltonian [tends to it in the Ax — 
limit] and gauge-invariant [accounts for the gauge trans- 
formation A — > A-|- Vx(i") inducing wave function phase 
change 'l'(r) — >■ exp (— ^x(r)) *(r)]. The kinetic energy 
operator— acting on a wave function defined on a mesh 



^^^'"^ = 2m* Ax^ ^"^^^'^ ~ ^^^^''--^ ~ 

= E■^^.,u■ (3) 

We find the energy E by solution of the boundary prob- 
lem in the incoming lead (see the next section), then Eq. 
is solved as a system of linear equations. 



C. Boundary conditions 

The confinement potential in both the input and the 
output channels depends only on the x coordinate. The 
chosen gauge allows for separation of the x and y coor- 
dinates in the Hamiltonian eigenfunctions 



^'(x,y) = exp(ifc2/)'0^(x). 



(4) 



with the wave vector k. In the absence of the magnetic 
field the nth Hamiltonian eigenstate of a w = 68 nm wide 
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channel has the energy £"„ = n^-K^ /2m*w'^ — \.21n^ 
meV. We consider the transport limited to the lowest 
n = 1 subband and skip the subscript n in the following. 
Only the electrons with wave vector exceeding k ~ 0.08 
nm^^ have enough energy to be scattered to higher sub- 
bands. We restrict our discussion to lower values of k. 
According to the Landauer-Biittiker approaches in the 
single subband transport the conductance G is simply 

2 2 

proportional to the transfer probability G = 

We first determine the boundary conditions for the in- 
coming lead. We assume v = 1 (lowest row of the mesh 
in the computational box) and plug 

*M,>^±i cxp(±ifcAx)V'^ (5) 
into Eq. ([3]) to obtain one-dimensional eigenequation 



2m*Aa;2 



cos(A:Aa:; + -BxAx) j 



(6) 



Eq. ([5]) provides the energy E that is used in the main 
equation ^ as well as eigenfunctions corresponding to 
incident (fc > 0) and reflected {k < 0) electrons that are 
used for setting the Dirichlet boundary condition for Eq. 
(U) at the bottom of the computational box 



(7) 



where the amplitudes of the incident Ck and reflected C-k 
wave functions are determined in a manner described in 
Section HTH 

Boundary condition ([7]) guarantees that the energy 
density H'^{x,y) /^{x,y) = E m the incoming lead and 
within the ring as found from Eq. ([3]) are equal. In 
order to match the energy density inside the ring the 
wave vectors in the output channels [see Eq. (|H)] in 
nonzero B must be different than k. Within the chan- 
nels the confinement potential is zero, therefore equal 
energy density means equal kinetic energy density. The 
kinetic energy operator is proportional to the square of 
the kinetic momentum = (p -I- eA)^ = Ul + Uy = 

(~*^^ ~'~ ■ Since the input and the out- 

put channels have the same width, the energy density is 
matched for the wave vectors in the left ki and right k^ 
output leads related to the wave vector of the incoming 
lead k as ki = k — ^xi, and kr = k — ^Xr, where xi 
and Xr arc positions of the axes of the left and right out- 
put leads {xr = —xi = 200 nm). Accordingly, for the 
boundary condition at the top of the computational box 
we use 



cxp(ifc'A2;), 



(8) 



where k' = fc; for a; < and k' = kr for a; > 0. This 
condition is introduced into Eq. for the top end of 
the computational box (v = 241). 



On the left and right edges of the computational box 
we introduce an infinite potential barrier which amounts 
in putting ^'^-i,^ or ^'^+1,,/ = in Eq. (jS]) for a mesh 
points at the left and right ends of the box, respectively. 
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FIG. 2: Probability density (lower panel) and probability den- 
sity current (upper panel) across the incoming lead for the 
incident ^''°'(a-) and backscattered \['^''''' (a;) electron eigen- 
functions for Ifcl = 0.05 nm"^ at B = 1 T. 




-50-40-30-20-10 10 20 30 40 50 
X rnml 



FIG. 3: The probability density in the incoming lead 
Ck'^'^ix) + C-k'^ ''{x) in the initial guess Ck = c-k and in 
the subsequent iterations of the self-consistent procedure (see 
text). Parameters are same as in Fig. [2] 



D. Backscattering probability 

The vertical component of the probability density cur- 
rent in the incoming lead 



j{x) = Q(<I»*-_) H Ay, 

TO* oy m* 



(9) 



is a superposition j{x) = j^{x) + 7 ^{x) of the incident 
current 

j\x) = A|cfc|2|^fc(a;)|2(^fc + ^Bx) (10) 
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and the backscattered one 



r'ix) - \c_k\'\i>-'ix)\'{-hk + eBx). (11) 

TO* 

Figure [2] shows the probabihty density and probabihty 
density current across the incoming lead for the incident 
(fc) and reflected (— fc) waves for k = 0.05 nm~^ and 
S = 1 T. Probability densities arc shifted from the axis 
of the lead to the left with respect to the direction of 
the current flow in consistence with the Lorentz force 
orientation. The basckattering probability is evaluated 
as a ratio of the current fluxes integrated across the input 
channel 



R = 



J dxj ''{x) 
J dxi'^{x) 



(12) 



For the axis of the incoming lead a; = 0, the solutions 
of the eigenequation (O with opposite k are related as 
il)^{x) = 'ip~^{—x) which implies that 1) the backscatter- 
ing probability is simply 



R 



C-k 



(13) 



and 2) both the incident and reflected wave function cor- 
respond to the same average value of H?. 

Solution of the system of equations (|o]) gives the wave 
function in the entire system. Now our task is to ex- 
tract c±k, i.e. the contributions of the incident and 
backscattered wave functions. For that purpose we con- 
sider two points in the incoming lead near the bottom 
of the computational box. We typically take two low- 
est points of the axis of the lead (/i = 120, v = 1) and 
(/i = 120, 1/ = 2), the results are not affected by a specific 
choice of these points. Wave function for v = 1 is given 
by Eq. ([7]) and for = 2 we have 



Ck^Jj^ cxp(ifcAx) + c-ki>n exp{-ikAx). (14) 



The eigcnfunctions -i/^, ijj^ ^ are determined from Eq. ([6|). 
Formulas ([7]) and (|14p form the system of equations for 
Ck and C-k- 



E. Self-consistence for the amplitudes of the 
incident and reflected wave functions 



For nonzero B the Hamiltonian ([6]) depends on the 
sign of the wave vector, and the eigenfunctions for ±fc 
are different. We need to assume some initial values for 
Ck and c^k to set the boundary condition ([7]) for the sys- 
tem of equations (jSj . Solution of Eq. Q gives the wave 
function in the entire computational box, including the 
incoming lead, of which Ck and c-k can be extracted. 
The procedure to determine Ck and c_fc is performed in a 
self-consistent iteration with Ck = c^k = ^ assumed as 
the initial guess. The iteration converges quite fast. Fig. 
[3] shows the charge density across the incoming lead for 



k = 0.05 nm"-'^ at i? = 1 T. The final result differs con- 
siderably from the initial guess but the results of the sec- 
ond iteration only slightly differs from the first one. For 
parameters applied in Fig. [3] the result for the backscat- 
tering probability R= |c_fep/|cfcp converges from 1 (for 
the initial guess) to 0.021. Naturally, the iteration affects 
the results in the entire computational box. 

Note, that by the initial guess Ck = C-k one assumes 
that the wave function in the incoming lead is symmetric 
with respect to its axis. For nonzero B this is the case 
only when backscattering probability reaches 100%. 



F. Transfer probabilities to the left and right 
output channels 

We need to separate the electron transfer probability 
to the left Ti and right Tr leads of the total transfer 
probability T = 1 — R. For that purpose we calculate the 
probability currents in the left and right output lead at 
the top of the computational box. Formula (jH]) with the 
boundary condition ([5]) gives 



ji{x) = —\^{x,y')\^ {hki + eBx)) (15) 

TO* 

for the left lead, and 



jrix) = —\^ix,y')\'{hkr + eBx)) 



(16) 



for the right one, where y' is the coordinate of the top of 
the computational box. We integrate the current fluxes 
on the left and right sides of the box J/ = /_i20Ax 
and J,. = Jq ^^'^^ dxjr{x). The transfer probability to the 
left and right channels is then calculated as Ti = Tj-^ 
and Tr ~T 



Jl+Jr- 



Jl+Jr 



G. Time dependent simulations 

For the interpretation of the results it is useful to con- 
sider also the solution of the time-dependent Schrodinger 
equation, ihd'^/dt = H'i'. For the initial condition we 
use a Gaussian wave function entirely localized in the 
input lead 



Afci/2 



^ix,y,t = 0) = J^^M^y- 



■-{y-Yf+Iqy 



(17) 



where Y lies far enough below the ring in the incoming 
lead. Probability density of the initial condition in the 
wave vector space is 



|5f(/fc)|2 = Cexp (-2(fc - g)VAfc2 



(18) 



The time-dependent calculations arc performed using the 
Crank-Nicolson scheme with a time step of 0.3 fs. We 
use the finite difference Hamiltonian with the same 
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mesh spacings as in the time independent calculation, 
but with radically enlarged computational box. For the 
time dependent simulation the computational box that 
we use covers as much as 12 yum of the input and the 
output channels. Computational box of this large size 
was necessary for setting the initial condition for a nearly 
monoencrgctic wave packet. 



0.5 
0.4 
0.3 
0.2 
0.1 



with 



II 



0.0 





VA., 


1 


1 1 




B=0 


\l 






, 1 , 




1^ 



0.02 



0.03 



0.04 



a) 



0.05 
k[1/nm] 



0.06 



0.07 



0.08 




0.0200 



b) 



0.0400 0.0600 
k[1/nm] 



0.0800 



FIG. 4: Transfer probabilities to the left T; and right Tr 
output channels and their sum T as functions of the incident 
wave vector /c for B = (a), and B = 0.8 T (b). 
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FIG. 5; Zoom of a fragment of Fig. [IJb) 



Simulation of the temperature effects for 
broadening of the transport window 



In order to estimate the effects of non-zero temperature 
for the transport we apply the linear response formula for 
the conductance^ 



G 



2e2 



(19) 



T= lTiE)[-^]dE, 



(20) 



and the Fermi function / = (^Q{E-EF)/kBT _|_ 2) ^ ^ where 
T stands here for the temperature. Formula (PH)) ac- 
counts for averaging the transfer probability obtained in 
the Hamiltonian eigenstates within the transport window 
that is opened near the Fermi level by thermal excita- 
tions. In the integral over the energy the wave vector k 
corresponding to a given E is found from the eigenequa- 
tion (il). 



III. RESULTS AND DISCUSSION 

A. results for K 

Calculated transfer probabilities as functions of the 
wave vector are presented in Fig. 21 For B = one ob- 
tains Ti{k) = Tr{k) due to the symmetry of the structure 
[Fig. |4|a)]. Non-zero magnetic field introduces asymme- 
try in the transfer probabilities [Fig. lU^b)]. Generally, at 
B > Q one observes that Ti is enhanced at the expense of 
Tr which is consistent with the orientation of the Lorentz 
force. Nevertheless, for discrete values of k sharp dips of 
T appear at higher B [Fig. IH^b)]. The dips of T coincide 
with the minima of T/ and peaks of T^. A zoom of one 
of T dips is shown in Fig. [51 The amplitude of the wave 
function and probability current distribution for k near 
the dip are displayed in Fig. [6l For k ~ 0.057 nm~^ 
the electron is directed to the left arm of the ring and 
then to the left output channel as previously described 
by the time dependent calculations J^"— For k = 0.0577 
nm""'^ the current forms vortices in the left arm and the 
actual electron transfer occurs through the right arm of 
the ring. For k = 0.0577 nm~^ the current goes through 
the right arm but the electron transfer to the right lead 
has still low probability. For k = 0.0579 nm^^ - at the 
center of the Ti dip (T^ peak) - the current forms a giant 
counterclockwise vortex around the entire ring. A min- 
imum of the wave function amplitude is formed at the 
center of the entrance to the left output channel - similar 
to the one observed for k = 0.0577 nm~^ at the right 
output channel. For larger k the current restarts to flow 
through the left arm as guided by the classical magnetic 
forces. 

The magnetic forces influence the distribution of the 
charge density within the ring. In Fig. [Hlwe observe a dis- 
tinct shifts of the wave function amplitude with respect 
to the axes of the channels correlated with the direction 
of the current and consistent with the orientation of the 
Lorentz force. For k = 0.057 nm~^ the wave function 
is distinctly shifted to the left edge of the input and the 
output channels as well as to the external edge of the left 
arm of the ring. For k = 0.0577 nm~^, when the transfer 
of the current through the left arm is blocked, the wave 
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FIG. 6: The red contours show the absolute vahie of the wave function (the darker the shade of red - the larger |'I'|) and 
probability current field (arrows) for B = 0.8 T and several values of k indicated at the top of the figure. For the transfer 
probabilities see Fig. (5] 



function maxima between the input and left output lead 
are placed symmetrically between the internal and exter- 
nal edges of the ring. For the giant anticlockwise vortex 
found for k = 0.0579 nm~^ the wave function is pushed 
to the inner edge of the ring. 

In experiments the conductance is usually measured in 
function of the magnetic field. Fig. [Tl^a) shows the trans- 
fer probabilities as functions of B ior k = 0.0683 nm~^. 
At low B the maxima of T correspond to interlaced peaks 
of Ti and T^. At higher B the value of Ti increases on 
average and the peaks of Tr become very narrow. Pro- 
nounced dips of Ti arc formed at the positions of Tr max- 
ima. The peak/dip structure occurs periodically with the 
spacings of AB = 0.09 T which corresponds to the flux 
quantum threading the one-dimcnsional ring of an effec- 
tive radius 121 nm that well agrees with the geometry of 
the model structure (Fig. 1). 

In order to quantify the direction of the current flow 
within the ring we calculate the flux of the current at the 



horizontal cross section of the arms of the ring y = 240 
nm (see Fig. 1). The fluxes arc then normalized to obtain 
-\- = 1. In Fig. mb) we notice that for larger B 
outside T dips nearly all the current goes through the 
left arm of the ring. 

Fig. |5] corresponds to fc = 0.0667 nm~^ for which a 
maximum of T = 2Ti = 2Tr is found for i3 = [see Fig. 
lUJa)]. At low B the peaks of Ti and Tr appear very close 
to one another forming a wider T maxima. For higher 
B i) the maxima of Tr turn into narrow peaks which 
coincide with the dips of T; ii) outside the T dips the 
current flows up through the left arm of the ring while 
the current flux through the right arm is close to zero, as 
discussed above for k = 0.0683 nm~^. 

Enlarged fragments of Fig. [5] corresponding to two dips 
of Ti are shown in Fig. [9] The amplitude of the wave 
function and the probability density current for the B 
range of Fig. M)^) is illustrated in Fig. [TOl For B = 0.35 
T the transfer probability to the left lead is maximal, 
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FIG. 7: (a) Transfer probabilities to the left Ti and right 
Tr output channels and their sum T as functions of B for 
k — 0.0683 nm^^. (b) The flux of the current through the 
left and the right arms of the ring. 




FIG. 8; Same as Fig. [I]but for k = 0.0667 nm" 



while Tj. is minimal. The current goes nearly entirely by 
the left arm. Note the pronounced elongated minimum 
of the wave function at the exit to the right output ehan- 
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FIG. 9: Zoom of two fragments of Fig. [8] corresponding to Tr 
maxima. The solid lines show the transfer probabilities (left 
vertical axis), and the dotted ones the normalized flux of the 
probability density current through the left and right arms of 
the ring (right vertical axis). 



nel. This wave function would be effectively coupled to 
the second subband of the right channel, but the latter 
corresponds to a much higher energy, so the transfer to 
the right lead is blocked. For B = 0.3536 T a leakage of 
the current to the right lead is observed and Jr becomes 
equal to Ji. For B = 0.3576 T the vortices of the cur- 
rent appear in the left arm. The transfer of the current 
through the left arm is nearly blocked. Note the posi- 
tion of the sharp minimum of the wave function near the 
left output lead pointed by the blue line in Fig. [TOl As 
B grows from 0.35 T the minimum is shifted to the left 
and for B = 0.3625 T it is found at the center of the 
junction of the left output channel to the ring. For this 
value of the magnetic field Ti becomes minimal. Gener- 
ally in our simulations a minimum of the wave function 
amplitude at the center of the junction to the left output 
lead is found for all minima of Ti which become sharp at 
higher B. When the electron transfer to the left lead is 
blocked or hampered, the current goes to the right out- 
put channel leading to appearance of a maximum of T^. 
For B = 0.3625 T the current forms vortices in the left 
arm as well as between the output leads and the main 
electron transfer goes through the right arm to the right 
output lead. For B = 0.39 T the minimum of the wave 
function is shifted to the lower edge of the left junction 
and the current transfer through the left arm to the left 
output channel restarts. 

For Jr > Ji the direction of the current circulation 
is opposite to the one preferred by the Lorentz force. 
Intervals of B corresponding to this orientation of the 
current become narrow at higher field [see Figs. [7Jb) and 
[51[b)]. Also the magnetic field interval for which Tr > Ti 



FIG. 10: The contour plot shows the absolute value of the wave function (the darker the shade of red - the larger |^[) and 
probability current field (arrows) for k = 0.0667 /nm and values of the magnetic field of Fig. [Sjb). 



become narrower at higher field [see also Fig. IH]. 

In order to conclude this section we note, that at higher 
magnetic field the electron transfer goes predominantly 
through the left arm of the ring to the left output lead, as 
should be expected due to the orientation of the Lorentz 
force. For narrow intervals oik or B wave function inter- 
ference within the ring leads to formation of a wave func- 
tion minimum at the entrance to the left output channel 
which blocks the transfer to the left lead. The T; minima 
are associated with reversal of the current circulation and 
appearance of maxima which turn into sharp peaks at 
higher B. 



B. Wave packet simulation 

The results presented so far indicate that for some 
intervals of the magnetic field the current flows in the 
opposite direction to the one indicated by the Lorentz 
force. The results of the wave packet simulation for 
nearly definite values of the packet wave vector should 
provide the transfer probabilities close to the ones found 
for the Hamiltonian eigcnstatcs. However, by the Ehrcn- 



fcst theorem in the wave packet dynamics the average 
values of electron momentum and position follow classi- 
cal laws. Hence, for i? > a preferential injection of 
the packet into the left arm of the ring is should be ex- 
pected for any magnetic field, on the contrary to the 
anomalous current injection that is found for Hamilto- 
nian eigenstates for some values of B. In order to inspect 
this contradiction closer we performed wave packet sim- 
ulations, in which we assume Afc = 5.5 x 10~^ nm [see 
Eq. (|17p ]. This wave vector dispersion for the studied 
structure and kp = 0.0667 nm~^ corresponds roughly 
to the thermal widening of the transport window which 
occurs at 150 mK. The spatial spread of the initial wave 
function is then as large as 4 /xm and we localize the wave 
packet Y = —8 /xm below the ring in the initial condition 

[Eq.dnD]. 

Figure ITlT a.b) shows the parts of the wave packet in 
the leads and within the ring for B = 0.35 T and B = 
0.3625 T. In Fig. [TlT b) we notice an enhanced packet 
transfer to the right output lead in consistence with Fig. 
m^b). Fig. [T2] shows the snapshots of the wave function 
amplitude and the probability current distributions for 
B — 0.3625 T. When wave packet enters the ring more of 
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FIG. 11: Parts of the wave packet in the input and the output 
leads as well as within the ring for q = 0.0667 nm~^ in the 
magnetic field of B = 0.35 T (a) and 0.3625 T (b). In (c) 
the lines show the transfer probabilities as functions of the 
magnetic field obtained by a time-independent calculation in 
which the transfer probabilities are averaged over the Gaus- 
sian distribution corresponding to the wave packet [Eq. (|2H) ]. 
The dots show the results of the time dependent calculation: 
the parts of the wave packet that are found in the output 
leads at the end of simulation for both B considered in (a) 
and (b). 



the electron wave function goes into the left lead (t = 43.5 
ps and t = 49.3 ps). At t = 70 ps an elongated wave 
function minimum is found at the entrance to the left 
lead. The current flow to the right output lead is visibly 
enhanced. For t = 79 ps the part of the wave packet 
inside the ring is maximal and we find that both the wave 
function amplitude and the current distributions are very 
close to those found in the Hamiltonian cigenstate for 
k = 0.0667 nm [see Fig. ^iov B = 0.3525 T]. 

Summarizing, in the time dependent simulations with 
a nearly monoenergetic wave packet one first observes an 
asymmetric injection of the packet to the arms of the ring 
in accordance with the Lorentz force orientation. Next 
the interference conditions similar to the ones found in 
the Hamiltonian eigenstates are formed. For B = 0.3525 
T the interference blocks the electron transfer to the left 



lead. 

The presented results of the wave packet simulation 
were obtained for an extremely low value of A/c. The 
time dependent simulations are useful for observation of 
the enhanced electron transfer to the right lead only for 
relatively low values of B, before the Tr maxima turn into 
peaks as sharp as in Fig. [8l[a) for B ~ 0.8 T. The AA: 
applied here corresponds to roughly 1/7 of the length of 
horizontal axis of Fig. [SJ which largely exceeds the width 
of the Ti dip. 

In Fig. Illf c) we compared the transfer probabilities 
estimated by the wave packet simulation with the ones 
obtained by the time-independent approach after calcu- 
lating an averaged over the wave packet probability den- 
sity in k space, i.e., 



{T) = C J dkT{k) cxp {-2{k - q f/Ak^ 



(21) 



Fig. [TTT c) shows that the results of the wave packet simu- 
lations are consistent with the /c-vector averaged transfer 
probability as calculated for Hamiltonian eigenstates. 



C. Finite temperature effect 

At high magnetic field the interference conditions lead- 
ing to anomalous injection of the current to the right arm 
of the ring appear for narrow k intervals. The conduc- 
tance measurements are performed in finite temperatures 
of the order of 100 mK^^i^ for which a transport win- 
dow of a finite width is opened near the Fermi level. In 
order to study stability of these anomalous transport con- 
ditions in finite temperatures we performed calculations 
for averaged transfer probabilities according to Eq. (|20p . 

For the temperature r = 115 mK the weight function 
— ^ calculated for Fermi wave vector— kp = 0.0667 
nm~^ is nearly a Gaussian function of k centered ai kp 

and 

T. In the B ^ oo 
limit the energy tends to the lowest Landau level for any 
wave vector E{k) — > htOc/'i, hence the widening of the 
k window for a given thermal energy ki,T at higher B. 
For T = 350 mK (700 mK) the corresponding half widths 
are Ak = 1.4 x 10"^ nm^^ (3 x 10"^ nm"^) and Ak = 
1.6 X 10-3 nm-i (3.5 x 10"^ nm-^), for S = and 0.8 
T, respectively. 

Fig. [T3]shows the transfer probabilities and normalized 
current fluxes for Fermi wave vector fixed at kp = 0.0667 
nm"^ and three values of the temperature (results for K 
were given in Fig. [8|). In finite temperature the dips and 
peaks of the transfer probabilities are transformed into 
smooth cxtrcma of reduced amplitude which eventually 
disappear at high magnetic field. The attenuation of the 
Aharonov-Bohm oscillations of the transfer probabilities 
for non-zero temperatures at higher B is in agreement 
with the results of previous wave packet simulationsii"— , 
in which the averaging of the transfer probabilities with 
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FIG. 12: Snapshots of the time-dependent simulation for the average wave vector q — 0.0667 nm^^ and Afc — 5.5 x 10^"' nm^^ 
[see Eq. p7|) ] for chosen moments in time. The contour plots show the amplitude of the wave function and the arrows - the 
current distribution. The color scale for the amplitude is the same for all the plots. The scale for the current vectors is different 
in each plot. 



k are embedded in the initial condition. The attenua- 
tion was also observed in the experimental data of Rcf. 
[20j . Results of Fig. [T3] indicate that the oscillations of 
the direction of the current circulation around the ring, 
which determine the orientation of the generated mag- 
netic dipole moment, are more thermally stable than the 
oscillations of the transfer probabilities, which determine 
the conductance. 



D. Ring with a perturbed potential 

The experimental results (Fig. 1 of Rcf. 2^) indi- 
cate a significant anisotropy of the potential landscape 
within the ring since already at i? = the conductance 
of one of the output leads largely exceeds the other. The 
appearance of peaks of Tr at high B > Q that we dis- 
cussed above were associated with specific interference 
conditions for which the electron wave function at the 
junction to the left output lead possessed a minimum at 
the axis of the lead (see Fig. 10 for B = 0.3625 T for in- 
stance) . A question which seems natural is whether such 
interference conditions are still possible for a quantum 



ring containing a potential defect. 

In order to answer this question we considered a 
perturbation introduced by Gaussian potential Vd = 
M^exp(— [(a; — + iv ^ yc)^]/^^], centered in point 
Xc = —104.8 nm j/c = 179.5 nm in the left arm just in 
between the input and left output leads. The size of the 
defect is assumed Rd = 30 nm. The results for the trans- 
fer probabilities and current fluxes are displayed in Fig. 

M 

For W = —5 meV the impurity introduces a potential 
cavity which mainly shifts the phase of the wave function 
passing through the left arm [Fig. [HT e.h)]. We observe 
no pronounced effect for the qualitative features of the 
transfer probabilities at high magnetic field as compared 
to a clean ring W ~ case (cf. Fig. 8). 

A potential barrier that is introduced for > ham- 
pers the electron transfer through the left arm. For 
W > 4 meV the transfer probabilities to the left and 
right output lead become distinctly different near B = 
[Fig. [HT c.e)]. the amplitude of the Aharonov-Bohm os- 
cillation is significantly reduced and the peak / dip struc- 
tures disappear in the high field limit. 

Results of Fig. [T3]for W = i and 5 meV resemble the 
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FIG. 13: (a,c,e) Transfer probabilities to the left T; and right output lead Tr as well as their sum T averaged over the thermally 
widened transport window for for kp = 0.0667 nm^^. (b,d,f) Normalized current fluxes through the left and right arm of the 
ring. The results are presented for the temperatures r = 115 mK (a,b) and 350 mK (c,d) and 700 mK (e,f). 



measured conductance.— Near _B = the electron trans- 
fer goes mainly to the right lead. Ti exceeds Tr only for 
B > 0.5 T. Note, that also for B > 0.5 T the current 
flux through the right arm of the ring largely exceeds 
the one through the left arm [Fig. [Hfd.f)]. The domi- 
nant electron trajectory for this transport conditions was 
indicated in Ref. [loj using wave packet simulations. 

We conclude that the presence of a repulsive potential 
defect induces not only the asymmetry of the transfer at 
B = and a weak amplitude of the Aharonov-Bohm os- 
cillation but also the absence of Tr peaks at high B. For a 
strongly asymmetric potential the peaks of Tr disappear 
also in zero temperature. 



E. Ring of an increased channel width 

The above results were obtained for the width of the 
channel within the ring fitted to the width of terminals. 
For an increased width of the ring channel the electron 
coming of the lowest subband of the input lead may pos- 
sess enough energy to occupy locally - i.e. within the ring 
- the second subband. A local scattering to the second 



subband may influence the mechanism of the electron 
transfer through the system. In order to study this point 
we decreased the inner radius of the ring from 88 to 68 
nm (see the inset to Fig. [T^ . 



The fc-resolved transfer probabilities are plotted for 
B = 0.8 T in Fig. [TH For fc < 0.06 nm we flnd similar 
results to the ones presented above: the transfer goes to 
the left output lead for nearly each value k. For k > 0.06 
nm~^ the scattering to the second subband of the channel 
becomes allowed and one observes a non-regular depen- 
dence of the transfer probabilities, with Tr exceeding T/ 
on some intervals. For k = 0.05 nm~^ [Fig. [rSTb)] the 
transfer probabilities change with the magnetic field in 
the same manner as for the ring of smaller width. Very 
different results are obtained for k = 0.0667 nm ~^ [Fig. 
[TSTc)]. Not a sign of periodicity can be noticed. The re- 
sults seem chaotic with no clear signature of the Lorentz 
force effect. The results for k > 0.06 nm~^ resemble 
rather the transport through a chaotic cavity (quantum 
billiard^l) than through a quantum ring. 
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FIG. 14: The transfer probabilities and normalized current fluxes for a repulsive potential defect of height 3 meV (a,b), 4 meV 
(c,d) and 5 meV (e,f). Plots (g,h) correspond to an attractive defect of depth —5 meV. 



IV. SUMMARY AND CONCLUSIONS 

We have discussed the role of magnetic forces in sta- 
tionary electron flow through a three-terminal quantum 
ring as obtained for Hamiltonian eigenstates in a single 
subband transport regime. We have shown that in most 
cases at high magnetic field the transport seems governed 
by the magnetic forces: the entire current is injected into 
the left {B > 0) arm of the ring and then ejected to the 
left output lead, with the transfer probability that tends 
to 100% at high magnetic fields. Exception to this rule 
are found only for narrow windows of magnetic fields for 
which interference conditions within the ring lead to for- 
mation of wave functions which are weakly coupled to the 
left output channel. This form of interference is associ- 
ated with anticlockwise circulation of the current within 
the ring and with an appearance of narrow peaks of trans- 
fer probabilities to the right output lead. The anticlock- 
wise circulation is anomalous from the point of view of 



the direction of classical magnetic forces since the current 
is injected into the right and not the left arm of the ring. 
The sharp peaks of the transfer probability to the right 
output lead that are found for high B disappear in finite 
temperatures for which the Aharonov-Bohm oscillations 
of conductance arc eventually attenuated. Oscillations of 
the current circulation turn out to be more resistant to 
the thermal widening of the Fermi level than the transfer 
probabilities. We have demonstrated that the imbalance 
of the transfer probabilities at B = as well as the reduc- 
tion of the conductance oscillations that are introduced 
by a scattering center within the ring are associated with 
removal of the interference conditions leading to appear- 
ance of the peaks of Tr at high magnetic field. Wc have 
considered the ring with a width larger than the width 
of the channels. We demonstrated that the for wave vec- 
tors which allow for appearance of local scattering to an 
excited subband within the ring channel - the results for 
conductance become chaotic in function of B without a 



13 



1.0 



Q. 0.5 



a) 



b) 



J 



0.0 
0.02 



B=0.8 T 





clear signature of either the Lorcntz force effect or the 
Aharonov-Bohni osciUations. 



0.04 0.06 
k[1/nm] 



0.08 




0.4 0.6 
B[T] 





0.4 0.6 
B[T] 



0.8 



1.0 



FIG. 15: (a) Wave vector resolved transfer probabilities to the 
left and right output leads for the inner ring radius decreased 
from 88 nm (grey circle in the inset) to 68 nm (the black 
circle inside the grey one). Results for the ring channel width 
equal to the with of the lead channels were presented in Fig. 
Ill^b). Transfer probabilities as functions of B axe plotted for 
k = 0.05"^ and k = 0.0667 nm ^ in (b) and (c), respectively. 



Acknowledgements This work was supported by the 
AGH UST project 11.11.220.01 "Basic and applied re- 
search in nuclear and solid state physics" . Calcula- 
tions were performed in ACK-CYFRONET-AGH on the 
RackServer Zeus. 



^ Y. Gefen, Y. Imry, and M. Y. Azbel, Phys. Rev. Lett. 

52, 129 (1984); M. Biittiker, Y. Imry, and M. Y. Azbel, 

Phys. Rev. A 30, 1982 (1984); M. Biittiker, Y. Imry, R. 

Landauer, and S. Pinhas, Phys. Rev. B 31, 6207 (1985). 
2 R. A. Webb, S. Washburn, C. P. Umbach, and R. B. Lai- 

bowitz, Phys. Rev. Lett. 54, 2696 (1985). 
^ G. Timp, A. M. Chang, J. E. Cunningham, T. Y. Chang, 

P. Mankiewich, R. Behringer, and R. E. Howard, Phys. 

Rev. Lett. 58, 2814 (1987). 
* A. Fuhrer, S. Liischer, T. Ihn, T. Heinzel, K. Ensslin, W. 

Wegscheider, and M. Bichler, Nature (London) 413, 822 

(2001). 

^ S. Pedersen, A. E. Hansen, A. Kristensen, C. B. Sorensen, 
and P. E. Lindelof Phys. Rev. B 61, 5457 (2000). 
U. F. Keyser, C. Fiihner, S. Borck, R. J. Hang, M. Bichler, 



G. Abstreiter, and W. Wegscheider, Phys. Rev. Lett. 90, 
196601 (2003). 

^ W. G. van der Wiel, Yu. V. Nazarov, S. De Franceschi, 
T. Fujisawa, J. M. Elzerman, E. W. G. M. Huizeling, S. 
Tarucha, and L. P. Kouwenhoven, Phys. Rev. B, 67 033307 
(2003). 

* A. Miihle, W. Wegscheider, and R. J. Haug, Appl. Phys. 

Lett. 91, 133116 (2007). 
^ Y. Aharonov and D. Bohm, Phys. Rev. 115, 485 (1959). 
" F. Martins, B. Hackens, M.G. Pala, T. Ouisse, H. Sellier, 

X. Wallart, S. BoUaert, A. Gappy, J. Chevrier, V. Bayot, 

and S. Huant, Phys. Rev. Lett. 99, 136807 (2007); M. G. 

Pala, S. Baltazar, F. Martins, B. Hackens, H. Sellier, T. 

Ouisse, V. Bayot, S. Huant, Nanotechnology 20, 264021 

(2009). 



14 



T. Bergsten, T. Kobayashi, Y. Sekine, and J. Nitta, Phys. 
Rev. Lett. 97, 196803 (2006). 

R. Leturcq, L. Schmid, K. Ensslin, Y. Meir, D.C. DriscoU, 
and A. C. Gossard, Phys. Rev. Lett. 95, 126603 (2005). 
B. Grbic, R. Leturcq, T. Ihn, K. Ensslin, D. Renter, and 

A. D. Wieck, Phys. Rev. Lett. 99, 176803 (2007). 

" R. Leturcq, D. Sanchez, G. Gotz, T. Ihn, K. Ensslin, D. C. 
DriscoU, and A. C. Gossard, Phys. Rev. Lett. 96, 126801 
(2006). 

S. Gustavsson, R. Leturcq R, M. Studer, T. Ihn T, K. 
Ensslin, D.C. DriscoU, and A.C. Gossard A. C., Nano Lett., 
8, 2547 (2008). 

^® T. Ihn, M. Sigirst, K. EnssUn, W. Wegscheider and M. 
Reinwald, New J. Phys. 9, 111 (2007). 

B. Szafran and P.M. Peeters, Phys. Rev. B 72, 165301 
(2005); A. Chaves, G.A. Farias, P.M. Peeters, and B. 
Szafran, Phys. Rev. B 80, 125331 (2009); T. Chwiej and 
K. Kutorasinski, arXiv:1003.1861 

B. Szafran and P.M. Peeters, Europhys. Lett. 70, 810 
(2005). 

B. Szafran and M. Poniedzialek, Phys. Rev. B 80, 155334 
(2009). 

E. Strambini, V. Piazza, G. Biasiol, L. Sorba, and F. Bel- 
tram, Phys. Rev. B 79, 195443 (2009). 
P. Foldi, O. Kalman, M. G. Benedict, and P.M. Peeters, 
Phys. Rev. B 73, 155325 (2006). 

C. H. Wu and D. Ramamurthy, Phys. Rev. B 65, 075313 
(2002). 

P. Foldi, O. Kalman, M. G. Benedict, and P.M. Peeters, 



Nano Lett. 8, 2556 (2008); P. Foldi, O. Kalman, and P.M. 

Peeters, Phys. Rev. B 80, 125324 (2009). 

H. Pothier, P. Lafarge, C. Urbina, D. Esteve, and M.H. 

Devoret, Europhys. Lett. 17, 249 (1992). 

V. Chandrasekhar, R. A. Webb, M. J. Brady, M. B. 

Ketchen, W. J. Gallagher, and A. Kleinsasser, Phys. Rev. 

Lett. 67, 3578 (1991); D. Mailly, C. Chapelier, and A. 

Benoit, Phys. Rev. Lett. 70, 2020 (1993). 

N. A. J. M. Kleemans, I. M. A. Bominaar-Silkens, V. M. 

Fomin, V. N. Gladilin, D. Granados, A. G. Taboada, J. M. 

Garcia, P. Offermans, U. Zeitler, P. C. M. Christianen, J. 

C. Maan, J. T. Devreese, and P. M. Koenraad, Phys. Rev. 

Lett. 99, 146808 (2007). 

K. Wilson, Phys. Rev. D 10, 2445 (1974) 

M. Governale and C. UngarelU, Phys. Rev. B 58, 7816 

(1998). 

M. Biittiker, Phys. Rev. B 38, 9375 (1988). 

S. Datta Electronic Transport in Mesoscopic Systems, 

Cambridge University Press (1995). 

In the present calculation the relation between the wave 
vectors and energy, kp and Ep in particular are given by 
Eq. p. 

M. Ferrier, L. Angers, A. C.H. Rowe, S. Gueron, H. Bouch- 
iat, C. Texier, G. Montambaux, and D. Mailly, Phys. Rev. 
Lett. 93, 246804 (2004). 

S. Liischer, T. Heinzel, K. Ensslin, W. Wegscheider, and 
M. Bichler, Phys. Rev. Lett. 86, 2118 (2001). 



